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Abstract 

An extension of the Kramers-Kronig method for treatment of polarized infrared reflectance spec- 
tra from the face of low-symmetry crystals, where directions of principal dielectric axes depend 
on frequency, is proposed. It is shown, how to obtain the frequency dependencies of the complex 
reflectivity tensor components, using three reflectance spectra measured for different directions of 
linear polarization of the incident wave, when reflected wave is immediately sent to the detector. 
The problem is formulated in a form of the system of integral equations, and effective numerical 
technique is found for solving it. The question of the further recovery of the complex dielectric tensor 
on the base of the reflectivity tensor, is discussed. The case of the monoclinic crystals is considered 
in details. An example of the application of the extended Kramers-Kronig method to the reflectance 
spectra from the (ac)-face of the single-crystal bismuth oxide with the monoclinic lattice is given. 



The Kramers-Kronig method (KK) along with the method of dispersional analysis (DA) is the 
most popular technique for treatment of reflectance spectra. As is known, it consists using the 
Kramers-Kronig integral transformations for recovery of frequency dependence of the phase 9 of 
the complex reflectivity r = \/ r Rex.pi9 on the base of the real reflectivity R measured in wide 
enough frequency range. Applicability of such a procedure follows from the fact that the real part 
of the analytical function lnr = (1/2) In R + i9, depends only on measured function R(uj), while its 
imaginary part is unknown function 9{u). So the function 9{uS) may be found by integration: 
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where integral with singularity must be treated in the sense of principal value Q. 

The next step is calculation of optical characteristics of matter. The normal incidence of elec- 
tromagnetic wave to semi-infinite homogenious medium from vacuum is usually discussed. In this 
case the complex dielectric function e may be calculated by formula: 
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that followes from the known Fresnel formulas. 

An important advantage of the KK method (in comparison with DA) is that it allows to obtain 
the complex dielectric function without model calculations, provided that reflectance spectra are 
known in wide enough frequency range. 

However, in the above form, the KK method is valid only when dielectric function of crystal 
and reflectivity coefficient can be considered as scalars. It is true when direction of electro-magnetic 
wave polarization coincides with one of the principal axes of the dielectric tensor e. In this case 
eigenvalue of tensor e corresponding to particular axis, plays the role of scalar e. This condition may 
be satisfied in entire spectral interval only for crystals with symmetry no lower than orthorombic. In 
the case of low symmetry crystals (with monoclinic or triclinic syngony) the dependence of direction 
of two or even three axes from frequency takes place (so called dispersion of dielectric axes). In 
monoclinic crystals only one axis is fixed in respect to crystalographic system of coordinates; two 
other rotate in the perpendicular plane. In triclinic crystals all axes have alternating direction. It 
doesn't allow to direct polarization of the incident radiation so that the vector of electric field and 
one of the principal axes of the dielectric tensor would be parallel. Hence the KK method (and the 
method of DA also) needs to be generalized. 

Generalization of the method DA for phonon spectra in monoclinic crystals was first proposed 
by Belousov and Pavinich 0, who introduced for each IR-active phonon an additional angular 
parameter which characterize direction of its dipole momentum and expressed tensor e through 
An + 3 parameters, where n - number of phonon modes. For phonon parameters recovery several 
spectra, obtained for different polarizations of incident waves, were used . In [|j we applied this 
method in another form. It was shown that it is sufficient to measure three different spectra; every 
additional spectrum is their linear combination. We have measured spectra for three angles of 
polarization equidistant at 45°. Both in || and in || the reflectivity was calculated for the most 
popular situation when only one (the incident) wave is polarized. In this paper we offer method 
for the recovery of the complex reflectivity tensor starting from the same three spectra, as in the 
method DA in form ||]. This generalized KK method may be applied for any symmetry crystals and 
may be used in situations, when directions of principal axes of the dielectric tensor are not known 
in advance, as is in a case of low symmetry crystals. 

Singularity may be eliminated by substitution of In R(£) in numerator by ln_R(^) — \nR(uj). 
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Calculation of the reflectivity tensor 




On assumption of medium linear response, the normal reflection of electromagnetic wave from the 
crystal surface is wholely described by the complex reflectivity tensor: 

(3) 

that associates complex amplitudes of incident (i) and reflected (r) waves: E r = rEj (for agreement 
with the notation in Q and || , we now assume that the wave vector of the incident wave is directed 
along the axis y). The components of the complex amplitude may be phase-shifted that corresponds 
to the elliptical polarization of wave. We will further assume that the medium is nongyrotropic (it 
means, for example, that no magnetic field is applied) and there is no spatial dispersion. In this case 
reflectivity tensor like dielectric tensor is symmetrical: r xz = r zx || fjj. 
Having set the direction of the incident wave polarization by angle x'- 



sinx 



we will obtain for the measured reflectivity: 
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R(X)= i F i2 i if 12 = \ r xxCOSX + r xz smx\ + \r xz cosx + r zz smx\ • (5) 

The problem is to find spectral dependencies of six functions which determine the complex 
tensor f using three reflectance spectra R(0°), i?(45°) and R(90°). The structure of the right side 
of expression (||) doesn't allow to separate known and unknown functions by taking of logarithm 
as is possible in isotropic situation. Apparently, solving the problem of complex reflectivity tensor 
recovery by simple integration is impossible. Nevertheless, the problem may be formulated in terms 
of integral equations as is shown below. 

Presence of two parts in (|5|) result from the fact that the reflected wave is generally elliptically 
polarized. If another polarizer (analyzer) were placed on the way of the wave, one of adding would 
be eliminated. But the scheme with only one polarizer is more preferable, because of the following 
reasons. Firstly, analyzer decreases intensity of light that falls onto detector and diminishes the 
signal-to-noise ratio. Secondly, the scheme discussed is the most typical for measurements of IR- 
spectra and often realized in commercial IR-spectrometers. At last, for the proposed variant of the 
method KK and for the method DA on monoclinic crystals || the same input spectra are used. 

Three equations for unknown functions can be obtained if one associate them with the measured 
reflectance spectra by formula (|5|): 

R(0°) = \r xx \ 2 + \r xz \ 2 , (6) 
2i?(45°) = \r xx + r xz \ 2 + \r xz + r zz \ 2 , (7) 
#(90°) = \r xz \ 2 + \r zz \ 2 . (8) 

As was mentioned, additional equations for other polarization directions won't add any informa- 
tion, because every new spectrum must be a combination of first three ones. This statement follows 
from type of dependence of reflectivity from angle x in accordance with (]5|). 

Three more equations may be obtained from the KK equations for each component of the tensor 

r: 
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where all integrals are implied in the sense of principal value. Shift by one in the last equation is 
introduced because r xz , unlike r xx and r zz , as function of the complex frequency may be equal to 
zero in the upper semiplane. It disturbs analiticity of function lnr xz and results in inapplicability 
of KK transformations. In contrast, function ln(l + r xz ) is analytical in upper semiplane, because 
\r xz \ < 1. 

Equations (H)-© along with (p|)-(|lT|) constitute complete system of equation, that may be solved 
numerically. 



Calculation of the dielectric tensor 

Having obtained frequency dependencies of all components of reflectance tensor, we must determine 
optical characteristics of crystal, particularly, the complex dielectric tensor e. 

Let us make general speculations, valid for all types of crystals, without any suggestion about 
orientation of dielectric axes. At first, let's remember that tensor e in the center of Brillouin zone 
(at k — ► 0) depends on direction of k. It is governed by fact, that if wave polarization vector P has 
longtudinal (relatively to k) component Pm, then a macroscopic electric field appears inside crystal 
Emacro = — 4vrP|| = — 47rk(kP)/|k| 2 , which depend non-analytically from k at k — > [||. When 
the incident wave is normal to the surface of crystal the wave vector of induced polarization waves 
(polaritons) inside crystal will also be normal to the surface, because in accordance to boundary 
conditions the tangential component of k must be constant. That's why all obtained values of 
tensor components will correspond to this direction of k. 

The knowledge of the reflectivity tensor r with dimensions 2x2, for one of the crystal faces is 
not enough to recover 3-dimensional tensor e. Here we can obtain only partial information about e. 
There is frequently used definition of the transverse reflectance tensor ej_ , that correlate the vector D 
with the transverse component of vector E: D = ej_Ej_. Because vector D is always perpendicular 
to the surface, the tensor e± is two-dimensional. In accordance to boundary conditions, the vector 
Ej_ must be constant on both sides of surface. Hence, it is the tensor that define reflectivity 
tensor r. In g were derived formulas that show the relation among components of dielectric tensor 
and that of reflectance tensor. In the most compact form this relation may be expressed by formula: 



(i-f)(i + f)- a , (12) 



that is similar to (|2|) (1 is the unit tensor). Thus, in general case the KK method allows to recover 
tensor e±. The further analysis depends on the crystal symmetry and particular model, that define 
features of dielectric function. Relations between tensors e and are described in details in pj. 

Let us consider in details crystal with monoclinic lattice. In this case the IR-radiation reflectance 
from the monoclinic (ac)-plane doesn't reduce to the scalar case, because directions of two dielectric 
axes in this plane are frequency-dependent. Symmetry allows to extract from the tensor e the 
two-dimensional minor c ac , that determine dielectric properties in this plane. It coinsides with the 
transverse tensor e± and therefore is defined by formula (12). Only phonon modes B n that have 
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dipole momentum in (ac)-plane contribute to e ac . In the model of polarized Lorenz oscillators this 
tensor equals: 

s (, A - i°° 4. ( cos2 6 3 cos B i sin °i \ SjUj2 i n n 

e ac (u) — e ac + a ■ a -2a ]~ 2 2 : ' 

^-r y cos C/j sin Uj sin Pj y — w z — Z7ju; 

where 5j - oscillator strength, uj - transverse frequency, jj - damping coefficient, 9j - angle between 
direction of dipole momentum and the axis x, - high frequency dielectric tensor. Summation is 
performed over all B^-modes. 

It is convenient to use e ac = Sp e ac = e xx + e zz as uniform scalar characteristics for dielectric 
permeabilty in the (ac)-plane. In the model ( |i~3| ) it doesn't depend on directions of oscillator dipole 
momenta, in contrast to e xx and e zz apart: 

e ac = S P 6gg + ^ 2 Sj f . . (14) 



Function (14) correspond to similar dependency of dielectric function from phonon mode char- 
acteristics in scalar case. Also, similar to a well-known rule for scalar case, transverse frequencies of 
modes B u may be determined as positions of maxima on frequency dependency of Im e ac . This rule, 
of course, is valid if modes doesn't strongly overlap (interval between two neighbour frequencies ujj 
greater than their widths) . Under this assumption the polarization angles 9j (defined with ambiguity 
to 7r) may be approximately obtained from equation: 

cos * ; m e -;^;-; m , ( i 5) 

Im e xx (ujj) + Im e zz {uij) 



which might be drawn out if from imaginary part of sum ( |13| ) the contribution of only j-th mode 
(dominating at u = Uj) is considered. The choice of right root is determined by sign of off-diagonal 
element of Im e xz (uj ) . 

Features of numerical procedure 

Having substituted: 

ae ia = r xx , be ifi = r zz , ce^ = l + r xz , (16) 
A = R{0°), B = R{90°), C = R(45°)-(A + B)/2, (17) 
we'll find that system of equations changes to: 

A = a 2 + c 2 -2 cos 7 + 1, (18) 
£ = 6 2 + c 2 -2cos7 + l, (19) 

C = accos(a — 7) + 6ccos(/3 — 7) — a cos a — b cos (3, (20) 
2uj f°° lna(f) 

2oj f°° In &(£) 

IT Jo C 2 - ^ 2 

IT Jo £ — ^ 

with six unknown functions a, b, c, a, j3 7 and three known functions A, B C. 



2uj f°° lna(f) , , 

a w = / -^4^, 21 

2u f°° In 6(f) , . 

ft") = -—l. aZT^' (22) 
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Let's vary unknown functions so that they satisfy system of equations. Inside measured frequency 
interval we choose quite dense grid of base points u\ < UJ2 < ■ ■ ■ < w„. We approximately specify 
functions a, b c by they values a, = a(oJi), bi = b(u>i) and c, = c{uii), which assume as unknown 
variables. Then we calculate values of functions a, j3 7 in base points, approximately substituting 
in formulas (|2~i"|)-(|23"|) integration by summation: 



a i = ^ f a i jhxaj, fa = ^V^ln^-, 7; = ermine,-, i = l,...,n, (24) 
j=i j=i j=i 

where coefficients Oij are expressed from u)\ , . . . , u n according to method of approximate integration 
(type of interpolation function, method of extrapolation outside interval). For example, if one 
substitute a partly-constant function instead of real function f(oj). 



/» = 

the following values of coefficients o"y will be obtained: 



uj < (u\ + u; 2 )/2 

(uj-t + Uj)/2 < uj < (ujj + uj j+1 )/2, j = 2, 

UJ > (U n -i + U> n )/2 , 



, n — 1 



(25) 
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2u)i + ujj + ujj + i 

In ^ ~ Uj 

2wj + cjj + Wj+i 
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1 

2, 
n 



(26) 



It corresponds to extrapolation of low- and high-frequency spectra by constant function, that is 
acceptable for dielectrics. This method was used to treating model and real spectra presented in 
this paper. 

In each point there are three measured quantities Ai = A(u>i), Bi = B{u)i) and Cj = C(ui). After 
definition of set of values ctj , bi and q and calculating with help of (^4[) ai , Pi and 7$ , the following 
functional can be composed: 



1 

3n 



E 



=1 L 



+ 



Bi-Bj 
5Bi 



+ 



a 



5a 



(27) 



where trial quantities A[, B[ and C- are expressed through o«, bi, Ci, ai, Pi, 7, in accordance with 
(P^|)-(p0|), and 5Ai, 5B{, 5Ci are meansquare noise deviations of Ai, Bi, d. 

Let's vary 3re quantities en, bi and Cj to find a minimum of functional ( ]27| ) as function of these 
quantities. As initial approximation it is reasonable to take a, = yA~i, bi = yBl, a = 1, what 
is exact decision in a case when the reflectivity tensor (||) is diagonal. Set of values that makes 
functional approximateli equal one may be considered as approximate decision. 

For functional minimization we have used the Marquardt-Levenberg nonlinear optimization 
method described in Partial derivatives of functional by variables aj, bi and Cj were calcu- 
lated in analytical form. It is advisable to vary the density of base points according to density 
of spectrum "features" for diminishing n. In our calculations maximum number of base points n 
was 450. Despite great number of variables, minimization process converged quickly enough (5-7 
iterations). 

Certainly, it is not so obvious, why must minimization process always lead to true result, because 
functional F in multi-dimensional space may have local minima that do not correspond to required 
solution. Besides, the noise in input spectra may cause instability of minimization process. However 
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our experience shows that minimum which corresponds to right decision, is usually so deep that 
minimization procedure converges just to it. Experimental noise (in reasonable limits) doesn't affect 
the final result very much. Nevertheless, the convergence of method and its stability to experimental 
noises require special attention. 

Application examples 

Before application of the discussed method KK to real experimental data, we tested it on a set of 
model spectra obtained with a help of the dispersion formulas for monoclinic lattices ||, taking 
different sets of parameters of B u phonon modes (frequencies, intensities, widths and polarization 
directions). To examine stability of the KK method to noise, additional noise with amplitude 5% 
was applied to initial spectra. Calculated by the KK method spectra for all components of the 
reflectivity tensor and the dielectric tensor were in a good agreement with that obtained directly by 
DA formulas. On the figure 1 results of one such test are showed. 

Application of the discussed method to real reflectance spectra of single crystal a-Bi203, that 
has monoclinic lattice is shown as an illustration. On the figure 2(a) three reflectance spectra from 
the face (ac) are presented. They were used as input data for the generalized KK method. Details 
of experiment and results of data treatment by the DA method are described in ||. 

Calculated by the KK method imaginary part of dielectric function e ac of bismuth oxide is 
presented on the figure 2(b) along with its frequency dependence calculated by the DA method. In 
general, a good match between these data may be noted. Several deviations, however, demonstrate 
that common description of dielectric function by Lorentz oscillators, isn't quite adeqiate. Inasmuch 
as discussion of optic properties of bismuth oxide isn't topic of this paper, we restrict ourselves by 
this notice. 

Authors are thankful to E.L.Kosarev, A.I.Kleev and E.R.Podolyak for discussion of problems 
concerning solving of system of integral equations. 
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Figure captures 

Fig.l. Results of the generalized KK method application to the model problem about reflection 
of IR-radiation from the (ac)-face of monoclinic crystal with two polarized in the (ac)-plane Lorentz 
oscillators, with characteristics: uj\ = 200 , S\ = 0.56, 71 = 10 , 9\ = 30°; u>2 = 450 -1 , 
S 2 = 0.79, 72 = 20 - 1 , 6 2 = 60°. Components e™: e xx = e zz = 5.0, e xz = 0.0. 

(a) Reflectance spectra at different polarization angles x, obtained according to formulas of 
dispersional analysis with added artificial noise 5%. 1-0° (shifted at 1), 2 - 45° (shifted at 0.5), 3 - 
90°. 

(b) Components of reflectance tensor r, obtained as a results of treatment spectra shown on 
figure la with help of KK method (dots) and directly with help of the DA formulas DA without 
additional noise (solid line). 1 - Im r xx , 2 - Re r xx , 3 - Im r zz (shifted at -1), 4 - Re r zz (shifted at 
-1), 5 - Im r xz (shifted at -2), 6 - Re r xz (shifted at -2.5). 

(c) The imaginary part of the dielectric function e ac , calculated on the base of reflectivity tensor 
(figure 2a). Dots - reflectivity tensor was recovered using the KK method. Solid line - "true" 
reflectivity tensor calculated directly by the DA formulas. 

Fig.2. 

(a) IR reflectance spectra at T=300 K from the (ac) face of monoclinic crystal a-Bi2 03 at 
different polarization angles x- 1 - 0° (shifted at 2), 2 - 45° (shifted at 1), 3 - 90°. 

(b) Frequency dependence Im e ac , recovered by two methods from initial spectra drawn at figure 
2a. Dots - function values obtained with a help of the generalized KK method. Solid line - results 
obtained with a help of the DA method || . 
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